Measuring deviations from a cosmological constant: 
a field-space parameterization 
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Most parameterizations of the dark energy equation of state do not reflect realistic underlying 
physical models. Here, we develop a relatively simple description of dark energy based on the dynam- 
ics of a scalar field which is exact in the limit that the equation of state approaches a cosmological 
constant, assuming some degree of smoothness of the potential. By introducing just two parameters 
defined in the configuration space of the field we are able to reproduce a wide class of quintessence 
models. We examine the observational constraints on these models as compared to linear evolution 
models, and show how priors in the field space translate into priors on observational parameters. 



Studies of dark energy [lj are hindered by the lack of 
specific alternatives to the cosmological constant model. 
Once we allow that the dark energy density might vary 
in time, how it evolves (most often parameterised by the 
equation of state w(z)) appears to be virtually a free 
function in most scenarios, making the models difficult 
to constrain. As a result, the dynamical dark energy 
evolution is often quantified instead by somewhat arbi- 
trary phenomenological models, such as assuming that 
w(z) evolves linearly in redshift or the scale factor. How- 
ever, such models generally have little motivation from 
fundamental physics, and can often act in ways (such as 
crossing the w = — 1 'phantom divide,' thereby violat- 
ing the weak energy condition) which can be difficult to 
explain from an underlying physics model. 

One commonly invoked physical model for dark energy 
is quintessence, where a scalar field <j> rolling down its 
potential comes to dominate the energy density of the 
Universe 0, In such models, the evolution of the 
equation of state is determined by the potential V{4>). 
It has been shown that any evolution of the equation of 
state could result from a suitably chosen potential Q]; 
however, it does not necessarily follow that all possible 
dark energy histories are equally likely. Here we examine 
the question of what priors on the functional space of 
w{z) result from the scalar field model, and develop a 
two parameter description which captures the behaviour 
of quintessence when its energy density is slowly evolving. 

If we take dark energy to be effectively described by a 
scalar field <fi, it seems reasonable to expect some level of 
smoothness in the configuration space of the field. For 
instance, we might expect that relevant variations of the 
Lagrangian to occur only over field-space distances of or- 
der one in some "smoothness scale" M. In the absence of 
univocal theoretical indications, there appears no strong 
motivation for choosing any particular scale for M; how- 
ever, the observational requirements of quintessence nat- 
urally associate it with the Planck mass. The general 



requirements for the potential V of a minimally coupled 
quintessence field are that it be of the size of the present 
energy density, V(<j) ) ~ MpH% (where M P = (87rG)~ 1/2 
is the reduced Planck mass) and its second derivative to 
be as small as the Hubble scale, V"(0o) S Hq-i so that 
the field is light enough to still be evolving today. We can 
introduce the smoothness scale M by writing the poten- 
tial as V{(j>) = MpH§f(cj)/M) and by assuming the "well 
behaved" function / and its derivatives to be at most of 
order one. The condition V"(<f>o) < Hq is then naturally 
satisfied if M >M P . 

The observational status of the dark energy equation 
of state provides the other basic motivation for this work. 
Recently, w(z) has been confined (g, LD, [8|, l2( to a region of 
the parameter space relatively close to the cosmological 
constant behaviour w = — 1. For a minimally coupled 
scalar field this corresponds to a slow- roll regime <fi 2 /2 <C 
V which, in turn, implies that the region A<f> probed by 
the field at late times is in fact small in Planck units. 

For these reasons, it is interesting to find a field-space 
parameterization of dark energy. In slow-roll inflation it 
is customary to assume that the dynamics is independent 
of the initial conditions and one defines field-space quan- 
tities (the slow-roll parameters) which capture the sub- 
stantial features of the cosmic evolution. The presence of 
a significant dark matter component makes things more 
difficult in the case of quintessence. Working in units 
2, the basic equations are 



where Mp 



4> + 3H4> + V'(<f>)=0, 
6H 2 = p m + p^, 



(1) 
(2) 



where p^ — 4> 2 /2 + V(<f>) and p m are the scalar and non- 
rclativistic matter energy densities respectively. The size 
of the acceleration of the field relative to the damping 
term, j3 = <p/3H<fc, is no longer systematically negligible. 
Typically, in thawing models discussed below, while the 
universe remains matter dominated (3 — 1/2, which then 
begins to decline towards zero when the quintessence field 
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begins to dominate (see also [10( about this) . By Eq. |T]) , 
the deviation from w = — 1 is given by 
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(3) 



Instead of using the usual slow-roll parameter e = 
(V'/V) 2 , we can try to reproduce the quintessence dy- 
namics by formally considering /3 as a function of <fi and 
making an approximation for the function 



V 



The following equations are easily found: 
1+w = k 2 0)^(1 - w) 2 /6 
= -/t(0)Q0(l - w), 



d\na 



(4) 



(5) 
(6) 



where = Pcf,/ (p m + P<p) and a = 1/(1 + z) is the scale 
factor. The above equations are as yet exact. When the 
dark energy density is nearly constant in time, Hw«l 
and Eq. ([SJ implies that k(</>) is also small. We can then 
consistently drop terms of order 0(1 + w) and we can 
substitute with its "zeroth order" approximation i.e. 
the fractional energy density corresponding to a cosmo- 
logical constant, 



f2 A (a) = 



(1 - fi0o)a 



(7) 



being the fractional energy density of dark energy to- 
day. (The differences in between fiA(a) and the true frac- 
tional density i~l^(a) are actually of order J a (l+w) da' /a' , 
but this should be no larger than 0(l+w) for the redshift 
range where the dark energy is dynamically interesting.) 
Equations and ([6]) can then be written as 



1 + w = 2n 2 ((j))n A (a)/3 
-2n((f>)n A (a). 
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(8) 
(9) 



The above equations, combined with an expansion for 
k(0), define our parameterization. There is still a one 
to one correspondence between the functions k(0) and 
w(a); given a particular k(4>) one can solve for w(a) and 
vice versa. However, the requirements of smoothness of 
the potential constrain the evolutions that are produced. 
While the full potential still has an infinite number of pa- 
rameters, the small field displacement when k(</>) is small 
ensures that only a narrow region of the field space is ac- 
tually probed by observations. The potential elsewhere 
is irrelevant; its associated parameters do not affect the 
likelihood and can be trivially integrated over in the evi- 
dence calculation, yielding an effective theory which has 
very few parameters. Note that the number of effective 



parameters of the model will depend on the values of the 
parameters themselves. For example, the possibility that 
the field sits at rest in a local minimum corresponds here 
to k(4>o) = 0. In this case, the model trivially reduces to 
a cosmological constant regardless of the values of k at 

<t> ^ 4>o- 

We found that a linear approximation for k(4>) = kq + 
Ki(<p — 0o) already captures very well the dynamics of 
models reasonably close to a cosmological constant. By 
inserting this linear form into © and integrating, we find 
4> — 4>o as a function of a. This gives 



k(q) = k (1 - + ^ 0o a 3 )- 2Kl/3 



K 



fi A (a) 



-i 2m/3 



(10) 

This and Eq. ([5]) specify uniquely w(a), which then can 
be integrated for the resulting evolution of the dark en- 
ergy density cx exp[T(a)], where 



1(a) 



:! I' % (1 + w(o!)) ~ 2 f % K 2 (a')n A (a>) 



= (K 2 (a) -«g)/2/si 



(11) 



Thus the dark energy density takes a very simple evo- 
lution in terms of Q\ (a). It is straightforward to check 
that our model is still well defined when k\ — > and k(o) 
is effectively constant. 

Quintessence models have been separated into two ba- 
sic classes 0]: those in which the field is effectively fixed 
at early times and only begins to evolve recently ( "thaw- 
ing" models) and those in which the field rolls quickly 
down a potential at early times, and only now is start- 
ing to slow down as the potential flattens out ( "freezing" 
models.) In the former, the equation of state begins at 
w = — 1 and rises at late times. In the latter, the equation 
of state begins at some higher value and drops towards 
w = —1 at late times. It is interesting to see how well 
our parameterization works for these models. 

In Figure [1] we plot the cosmological evolution of the 
dark energy equation of state for different values of Ko 
and K\. Our "smoothness requirement" is generally met 
for the potentials of the thawing models, and therefore 
our parameterization does very well at reproducing their 
behaviour. We find that the two parameters cover very 
well the entire range of the "thawing" quintessence mod- 
els. On the other hand, freezing models, while interest- 
ingly addressing the "coincidence problem" ll| , typically 
exhibit singular potentials and rather fine-tuned condi- 
tions for the field. Steep curvature and rapidly changing 
potentials at early times are in fact required in order to 
maintain a tracking [12| behaviour, and so are not well 
described by our parameterization. 

Our approach results in a two parameter (k , k\) fam- 
ily of w(a) histories which can be used to fit the data, 
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FIG. 1: The evolution of dark energy equation of state as 
a function of redshift, and in the (w, w') plane for different 
parameterizations. The solid (blue) lines show the results of 
an exact integration for a typical potential (quadratic in this 
case), while the dashed (red) lines show our parameterization. 
The approximation improves as the present equation of state 
approaches — 1. Also shown (dot-dashed) are similar curves 
for a linear parameterization, w(a) = wq + w a (l ~ a). We 
chose parameters to match w and w' at the present. 

much like the often-used linear parametrization, w(a) — 
wo+w a (l— a) I3,14|. If we use our parametrization to fit 
a fully dynamical thawing model, the values of (kq,k\) 
we obtain should reflect the true k(<P), modulo correc- 
tions of order 1 +w. 

To simplify comparison to the data and past literature 
[REi!]. we can transform from the (koj k i) parameters to 
w and w' = dw/dhia evaluated at the present: 



w{z 
dw 



0) 



2 Q^qKq 



(12) 
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-Kgf20 O [3 - ^0o(3 + 4/ti)] . 



A linear model with the same values of (w,w') at z = 



follows our model at low redshifts, but tends to diverge 
significantly by z = 1, as can be seen in Fig. 1. Note, 
matching at z — makes the differences between the 
models more apparent than if they were matched at a 
higher redshift; the model differences are generally small, 
particularly when translated into integrated observables 
like the luminosity distance. 

In Fig. [2] we compare directly the observational limits 
for our parameterization to those of the linear parameter- 
ization. We focus on data constraining the background 
evolution of the dark energy: the supernovae type la (SNe 
la) set of the SNLS survey Q as standard candles, the 
baryon oscillation (BAO) peak detected in the SDSS lu- 
minous red galaxies (LRGs) [l6| and the CMB shift pa- 
rameter from the WMAP 3-year data Q calculated by 
Wang and Mukherjee [13] as standard rulers. These di- 
rectly probe the cosmological expansion through the lu- 
minosity distance at low redshifts and the angular diam- 
eter distance at higher redshifts. The likelihood contours 
for the two parametrizations are very similar with the 
present data, reflecting their similarity in the low red- 
shift region where the observations are most sensitive. 

The likelihood curves in Fig. [2] ignore any prior in- 
formation about the models, effectively show the curves 
which would result from a flat prior in the (w, w') space. 
However, one might argue that it is more realistic to take 
priors in the space of potentials, e.g. assuming a uniform 
distribution in the derivatives of the potential or perhaps 
the log of the potential. In the limit of small /3, the lat- 
ter corresponds to a uniform prior on (kq,Ki). Thus it 
is interesting to investigate how this kind of prior would 
translate into the (w,w') space. 

This prior is straightforwardly given by the Jacobian 
of the transformation (|12p . which can be easily shown 
to be proportional to [fi<£o(l + w)]~ 3 / 2 . Thus, given a 
flat field space prior, the probability diverges as one ap- 
proaches w = — 1. While one might initially see such a 
prior as problematical, in fact this is not the case. The 
condition V"(<j>) < H$ effectively puts bounds on ki and 
restricts the interesting range of the (w,w') plane. Con- 
stant Ki translates into straight lines converging at the 
cosmological constant point, as shown in Figs. 2 and 3. 
Integrating the prior between such lines always results in 
finite answers. Also, while the prior will tend to pull the 
peak of the posterior towards w = —1, if the data clearly 
prefer a higher value, the divergence in the prior will be 
suppressed by the exponential falloff of the likelihood. 

Fig. [3] shows the prior and posterior for this parame- 
terization, assuming a flat prior on the (kq, k±) plane. In 
order to get confidence contours for the prior, we cut it 
off at w = —1/3, which corresponds to the border of the 
accelerating models and is also roughly where Kq = 1. 
Both the prior and the likelihood peak around w = — 1, 
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FIG. 2: We show contours of constant likelihood resulting 
from SN, CMB and baryon oscillation data. The upper figure 
assumes a linear model for the equation of state of the dark 
energy, while the bottom figure shows our scalar field moti- 
vated parameterization, where the space of models is reduced. 
The models have the same equation of state and derivative at 
z = 0, but evolve differently at higher redshift. For simplicity 
we have fixed Q^o = 0.74. 

making the posterior even more peaked. Our primary 
aim is to show how to translate priors from variables 
related to the theory to more observational parameters. 
The true nature of the prior on the (kq, k\) plane is still 
an open question and requires a better defined theory; 
the posterior shown in Fig. 3 should be only taken as 
representative of the kind of focussing that might occur. 

We have shown that as one approaches cosmological 
constant behaviour, the evolution of quintessence dark 
energy is more and more constrained by the requirement 
of a smooth potential, and the number of free parame- 
ters shrinks. This provides a unique template of models 
which can be used to search for small deviations from the 
cosmological constant behaviour. Our parameterization 
is exact in this limit, and so is more appropriate than a 



FIG. 3: The upper figure shows the prior on (w,w') space 
arising from a uniform prior in (/to, Ki) space, which diverges 
as (l+w)~ 3 / 2 . The points show how a uniform grid in (/to, Ki) 
is projected into the (w, w') space, while the lines show 1 and 
2cr confidence limits, assuming w < —1/3. The bottom panel 
shows the final posterior, where the prior is folded in with the 
likelihood from the data (from Fig. 2b.) 

linear model which can behave unphysically at high red- 
shift; other parameterizations, such as a subset of the 
'kink' models [l8| . would be an improvement to the lin- 
ear models, but they are still not exact and are harder to 
relate to the underlying theory. Our approach is similar 
in spirit to the flow equation approach of Huterer and 
Peiris [19(, which appears to find thawing models less 
likely. The reason for this discrepency is not obvious, 
but it likely lies in the initial conditions or in a different 
assumptions about the possible potentials. 

Since it is well known that the assumptions about w(z) 
in an analysis can dramatically affect one's conclusions 
about dark energy (2o| . it would be very interesting to 
repeat projections for future experiments using this more 
physical parameterization. It seems likely that, given the 
convergence towards w = — 1 at high redshift for our 
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parameterization, experiments at low and intermediate 
redshift will increase in importance. It would also be 
interesting to extend our approach, such as to models 
where the dark energy is coupled to dark matter 21 1 . 
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